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A five-dimensional Ricci-flat cosmological solution is studied by assuming that the in- 
duced 4D matter contains two components: the usual fluid for dark matter as well as 
baryons and a scalar field with an exponential potential for dark energy. With use of the 
phase-plane analysis it is shown that there exist two late-time attractors one of which 
corresponds to a universe dominated by the scalar field alone and the other is a scal- 
ing solution in which the energy density of the scalar field remains proportional to that 
of the dark matter. It is furthermore shown that for this 5D scaling solution the uni- 
verse expands with the same rate as in the 4D FRW models and not relies on which 4D 
hypersurface the universe is located in the 5D manifold. 
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1. Introduction 

Scalar fields play a central role in modern cosmology in driving inflation of the early 
universe and describing dark energy of the present epoch. 1,2 Observations predict 
that our universe contains roughly one-third of dark matter and baryons and two- 
thirds of dark energy. 3 Within the standard Friedmann- Robertson- Walker (FRW) 
models it was shown that there exist scaling solutions that are the unique late-time 
attractors whenever they exist. 4,5 ' 6 Here, in this paper, we wish to look for scaling 
solutions in higher-dimensional cosmological models. 

In Kaluza-Klein theories as well as in brane world scenarios, our 4D universe is 
believed to be embedded in a higher-dimensional manifold. One of these models is 
the Ricci-flat 5D cosmological solutions presented by Liu and Wesson. 7 This model 
is 5D Ricci-flat, implying that it is empty viewed from 5D. However, as is known 
from the induced matter theory, 8 ' 9 4D Einstein equations with matter could be 
recovered from 5D equations in apparent vacuum. This approach is guaranteed by 
Campbell's theorem that any solution of the Einstein equations in N-dimcnsions can 
be locally embedded in a Ricci-flat manifold of (N-fl)-dimensions. 10 In section II, 
we suppose the 4D induced matter be composed of a perfect fluid and a scalar field. 
In section III, we use phase-plane analysis to study the evolutions of the model. In 
section IV we study the scaling solution. Section V is a short discussion. 
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2. 5D Solution With An Induced Scalar Field 

An exact 5D cosmological solution was given firstly by Liu and Mashhoon in Ref. 
11 and restudied by Liu and Wesson in Ref. 7. This solution reads 

dr 2 



dS 2 = B 2 (t,y)dt 2 - A 2 (t,y) 
B = 



1 — kr 2 



r 2 (d6 2 



dy z 
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where A = (d/dt)A, /i = /j,(t) and v — v(t) are two integration functions of t , k is 
the 3-D curvature index (k = ±1, 0) and if is a constant. This solution satisfies the 
5D vacuum equation Rab = 0, with the three invariants being 

72K 2 

h = R = 0, 1% = R Rab — 0, I3 = R Rabcd 



.4 s 



So K determines the curvature of the 5D manifold. 
The 5D line element in contains the 4D one, 

T dr 2 

ds 2 = g aP dx a dx 13 = B 2 dt 2 - A 2 



\d9 2 +sin 2 



(2) 



(3) 



1 — kr 2 

Using this 4D metric we can calculate the 4D Einstein tensor by = ( - 4 - ) Rp — 5^ 



(4) 



R/2. Its non- vanishing components are 



(4)^0 — 3(p +k) 

- — a 2 — ■■ 



WGJ =W G 2 =< 4 ) G% = ^£ 



(p 2 + k) 
A 2 ■ 



(4) 



Generally speaking, the Einstein tensor in (0} can give a 4D effective or induced 
energy-momentum tensor via ^G^ = x?( 4 'TP with x 2 = 8ttG. 

In previous works, 7,12 this energy-momentum tensor was modelled by a perfect 
fluid with density p and pressure p, plus a variable cosmological term A. This A 
could be served to describe dark energy. In this paper, we let the energy-momentum 
tensor consist of two parts: 



p,v ' 1 fiV) 
^au ~ {Pni ~i~ Pmj'U'f/U'v Pra9^ivi 

Til = d u J)d u <t> ~ g^ih^d^d^ - V(4>% 



(5) 



where T™, represents a perfect fluid and T^ v represents a scalar field. Then, from 



the 4D conservation laws T 



(m) v 



and T 



(</>)„ 



0, we obtain 



Pm + ^J-iPm +Prn) = 0, 
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The Hubble parameter H and the deceleration parameter q should be defined in 
terms of the proper time as. 12 

AB A nA 
Using this we can recover the corresponding Friedmann equations as follows: 

H 2 + 7F = IT (Pm + P4>) 

2 

H = -\B{p m +p m + p,p + p$), 

where 

1 6 2 Id) 2 
P* = 2^+V(4>),P<f,= 552 (8) 

3. Phase-plane Analysis Of The Solution 

From (JTJ we see that in general we have B(t,y) ^ const, so t is not the proper 
time. This leads to the only difference between our results an d those of the 

standard FRW models. Similar as in Ref. 5, we define x and y as 



K(j) XVV 

y = -7^77- ( 9 ) 



V6BH V3H 

For an exponential potential V = Vq cxp(— \xd)) and a spatially flat universe [k = 
0), we find that the evolution equation for x and y are of the same form as in Ref. 

5, 

/3 3 
-y 2 + -x[2x 2 + 7m (l - x 2 - y 2 )}, 

V = - X \j\ x V + \v\^ + 7m(l - x 2 - y 2 )}, (10) 

where a prime denotes a derivative with respect to the logarithm of the scale factor, 
N = In A. Define 

n m = _|^,o = -^, (ii) 

then (|10|) gives fl m +0^ = 1 and JjJJ gives x 2 + y 2 — O^. Meanwhile, from |JSJ, the 
effective equation of state for the scalar field gives 

^ _ 2x 2 
P4> x + y z 

From the analysis of Ref. 5 we know that there are five fixed points (critical 
points) corresponding to x — and y = in the plane-autonomous system, in- 
cluding two stable nodes. 5 The two stable nodes represent two possible late-time 
attractor solutions and, therefore, are of particular physical interest. We list them 
in the following. 

(1). Late-time attractor solution dominated by the scalar field alone. 
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This solution corresponds to the stable node 

x = X/y/Q, y = V 1 - A 76, withX 2 < 6 (13) 
for which we have x 2 + y 2 — 1 , then we get 

n m = o,n^ = i,7 = a 2 /3. (u) 

So the universe is dominated by the scalar field alone and expands with a power-law. 
It can be shown that if, furthermore, A 2 < 2, the universe is accelerating. 
(2). Scaling solution. 
It corresponds to the stable node 



for which we have 



3 7m /3(2-7 m )7 m 2 . , 

2X' y V 2A 2 > mth * > 3 7m (15) 

A 2 - 3-y m 37 m 
"m = ^ ' il <i> = ~^2~>7<A = 7m- (16) 

So this represents a scaling solution where the dark energy density of the scalar 
field is proportional to that of the baryons and dark matter. 

4. 5D Scaling Solution 

It is known that in the standard spatially-flat FRW models the scaling solution (|15|) 
corresponds to a power-law expansion. Here let us check how the universe expands 
in the 5D model 0} . 

Firstly, by substituting the value of x of 1)15(1 into its definition in © and using 
(JZJ), we find that the scalar field 4> can be integrated, giving 

0= ^.lnA + d, (17) 

Thus we obtain the exponential potential as 

V = V Q exp(-A^) = C 2 V a A- 3 ^ (18) 

Meanwhile, substituting the value of y of H15|) into its definition in |[?JJ, and then 
using ||TBj, we get a constraint between A and /j, as 

fi 2 = C 2 A 2 - 3 ^. (19) 

Consider the 5D metric Q. The form Bdt = (A/fi)dt is invariant under an 
arbitrary coordinate transformation t — > tit). This would enable us to choose the 
coordinate t such that although t does not equal to the proper time r globally, 
it can tend to the proper time as t — > oo. Be aware that in the 4D standard 
FRW models the scaling solution 1(1 5JI corresponds to the expansion a(t) oc t 2 ^ 3lm ^ 
and H = 2/(3"/ m t). Thus we obtain the late-time approximation of the 5D scaling 
solution 1(15JI as: 

2 2 2 i n 2 , 

AMiJS.Bsl^s-tiS-U^MT-:. (20) 

37m ^4 37 m £ 
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For k = 0, the 5D solution JIJ reads 

v 2 + A 

Using J20J in (|2U, we find 

2 



A 2 = M V + 2^+ 2 . (21) 



37„ 



t^~\ (22) 



5. Discussion 

In this paper we have used the phase-plane analysis to study the stability of evolu- 
tion of a 5D cosmological model. From FCas. I|2f)fl - (|22fl we can draw three conclusions: 

(1) . The scaling solution (|15|) . as a late-time attractor of the evolution of the 
universe, is the same in both the 4D FRW models and the 5D induced matter 
theory. 

(2) . In the 5D model (|2"0"|) - (f2"2")l . the expansion rate A and H of (|2"0"|l seems not 
rely on the value of y. In other words, our 4D universe, no matter it is situated on 
which hypersurface of y in the 5D manifold, looks similar and expands with almost 
the same rate, in late times. 

(3) . The constant K in (|21|l represents the curvature of the 5D manifold (see 
Eq.©)- It is reasonable to believe that K = 0, ±1 may correspond to three different 
topologies for the 5D manifold However, from (|20|) - (|22ll . we find that K does 
not affect A and H so much in late-times. 

(4) . In Eq. (|2"Uj) . if ~f m — 1/3, we have B ps 1 and A w t 2 , and our model 
approaches to the Milne model which has no event or particle horizon, and so it is 
of particular interest. 

It was shown 7 ' 12 that in the early-time of the universe the 5D solution 
deviates from the standard FRW models greatly. For instance it may have a big 
bounce rather than a big bang and before the bounce the universe is contracting. 
Meanwhile, the value of the 5D curvature K is sensitive for having or not having 
a bounce. However, in the late time, we have B — > 1 and the coordinate time 
t gives back to the cosmic time r, so the bounce model approaches to the FRW 
model. As a late-time attractor, the scaling solution seems not show up noticeable 
deviations from the standard FRW models. Further study along this line is needed 
to distinguish the two theories. 
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